Resonant interactions of nonlinear water waves in a finite basin 
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We study exact four-wave resonances among gravity water waves in a square box with periodic 
boundary conditions. We show that these resonant quartets are linked with each other by shared 
Fourier modes in such a way that they form independent clusters. These clusters can be formed 
by two types of quartets: (1) angle-resonances which cannot directly cascade energy but which can 
redistribute it among the initially excited modes and (2) scale-resonances which are much more rare 
but which are the only ones that can transfer energy between different scales. We find such resonant 
quartets and their clusters numerically on the set of 1000 x 1000 modes, classify and quantify them 
and discuss consequences of the obtained cluster structure for the wavefield evolution. Finite box 
effects and associated resonant interaction among discrete wave modes appear to be important in 
most numerical and laboratory experiments on the deep water gravity waves, and our work is aimed 
at aiding the interpretation of the experimental and numerical data. 
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of energy in Fourier space and, therefore, they are analo- 
gous to the Kolmogorov energy-cascade spectrum in hy- 
drodynamic turbulence. WTT approach can also be used 
to study evolution of higher momenta of wave amplitudes 
and even their probability density function, and, there- 
fore, to examine conditions for deviation from Gaussian- 
ity and onset of intermittency @, 0, 0, HI- Significant 
effort has been done in the past to test WTT and its 
predictions numerically 0, @, 0, B as well as experi- 
mentally 13, 11, 12, HI- It has been noted however that 
in both the numerical and the laboratory experiments 
the domain boundaries typically play a very important 
role and the infinite-box limit assumed by WTT is not 
achieved. Indeed, as mentioned in 0], to overcome the 
wavenumber discreteness associated with the final box 
one needs at least 10000 x 10000 numerical resolution 
which is presently unavailable for this type of problems. 
On the other hand, as shown in (loT |. even in a 10m x 6m 
laboratory flume the finite-box effects are very strong. It 
is important to understand that WTT takes the infinite- 
box limit before the weak-nonlinearity limit, which phys- 
ically means that a lot of modes interact simultaneously 
if they are in quasi-resonance, i.e. satisfy the following 
conditions 



INTRODUCTION 



Weakly nonlinear systems of random waves are usu- 
ally studied in the framework of wave turbulence the- 
ory (WTT) (for introduction to WTT see book [K]). 
Besides weak nonlincarity and phase randomness, this 
statistical description is based on the infinite-box limit. 
This approach yields wave kinetic equations for the wave 
spectrum which have important stationary solutions, 
Kolmogorov-Zakharov (KZ) spectra [lj]. Importance of 
KZ spectra is in that they correspond to a constant flux 



|w(ki)±w(k 2 )±...±w(k 8 )| < n, 
ki ±k 2 ± ... ±k s = 



(1) 



with some resonance broadening > which is a mono- 
tonically increasing function of nonlincarity (mean wave 
amplitude). Here k and w(k) are wave vector and dis- 
persion function (frequency) which correspond to a gen- 
eral wave form ~ expi(kx — cut). WTT is supposed to 
work when the resonance broadening il is greater than 
the spacing <5 W between the adjacent wave modes 



n > {doj/dk)2Tr/L, 



(2) 
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where L is the box size. For some types of waves, for 
example for the capillary water waves, this condition is 
easy to satisfy, and therefore to achieve WTT regime 
[l4| However, this condition is often violated for some 
other types of waves, in particular for the surface gravity 
waves which will be the main object of this paper. This 
occurs in numerical simulations, due to limitations on the 
numerical resolution and in laboratory experiments, 
due to an insufficient basin size . In these cases, 
the Fourier space discreteness (which is due to a finite 
box size) leads to significant depletion of the number of 
wave resonances with respect to the infinite box limit. 
In turn, this results in a slowdown of the energy cascade 
through the fc-space with respective steepening of the 
wave spectra (lOl [Hj . In addition, the wavenumber grid 
will cause the wave spectra to be anisotropic in this case. 

What happens when the condition ^ is so badly vio- 
lated that only waves which are in exact resonance (i.e. 
fl = 0) can interact? In this case, the mechanism of 
the wave phase randomization based on many quasi- 
resonant waves interacting simultaneously will be absent 
and, therefore, one should expect less random and more 
coherent behavior. In [l6| it was shown that in many 
wave systems, resonantly interacting waves in finite do- 
mains are partitioned into small independent clusters in 
Fourier space, such that there cannot be an e nerg y flux 
between different clusters. In particular, in [171 ] some 
examples of wave systems were given in which no reso- 
nances exist (capillary water waves, u> — |k| 3 / 2 ), as well 
as systems with an infinite number of resonances (oceanic 
planetary waves, u> = Ikl" 1 ). Both of these examples are 
three- wave systems [i.e. s = 3 in ([!])]. In the present 
paper we will concentrate on finding resonances for the 
deep water gravity waves, which is a four-wave system. 

The problem of computing exact resonances in a con- 
fined laboratory experiment is highly non-trivial because 
the wavenumbers are integers and ([I]) is a system of Dio- 
phantinc equations on many integer variables in large 
powers. Computational time of solving this system by 
a simple enumeration of possibilities in this case grows 
exponentially with each variable and the size of spectral 
domain under consideration. A specially developed q- 
class method [IH, 0, [20| has allowed to accelerate the 
computation and find all the resonances among waves in 
large spectral domains, in a matter of minutes. In this 
paper we use g-class method to construct resonant wave 
clusters formed by 4-wave resonances among water grav- 
ity waves covered by the kinematic resonance conditions 
in the form: 

f|k 1 |V2 + |k 2 |V2 = 1^11/2 + ^11^ 
\ki+k 2 = k 3 + k 4 

with kj = (mi,ni) and integer |mj|,|n,-| < 1000. Our 
main aim is to understand how anisotropic resonance 



clusters influence the general dynamics of the complete 
wave field. 



CONSTRUCTION OF q-CLASSES 

We adopt the general definition of a g-class given in [l8| 
for the dispersion function u> = |k| 1//2 in the following 
way. Consider the set of algebraic numbers R = zb/c 1 / 4 . 
Any such number k has a unique representation 

k = 7<? 1/4 ,7 G Z, 

where q is a product 

q=p L 1 1 p e 2 2 ...p e n ™ 7 

while pi,...p n are all different primes and the powers 
ei,....e n <G N are all smaller than 4. Then the set of 
numbers from K having the same q is called a g-class 
Cl q . The number q is called a class index. For a number 
k = 7g 1 / 4 , 7 is called the weight of k. For instance, wave 
vector k = (160, 40) belongs to the g-class with q = 1700. 
Obviously, for any two numbers fci , k 2 belonging to the 
same g-class, all their linear combinations with integer 
coefficients belong to the same class q. 

It can be shown that Sys.Q have two general types of 
solutions: 

Type I: All 4 wave vectors belong to the same class Cl q , in 
which case first equation of the Sys.([3]) can be rewritten 
as 

71^ + 72^7 = 73^ + 74^7 (4) 

with integer 71 , 72 , 73 , 74 ■ 
and 

Type II: All 4 wave vectors belong to two different classes 
Cl qi ,Clq 2 ; in this case first equation of the Sys.([3|) can 
be rewritten as 

7i \fq\ + 72 \fq~2 = 7i \fk~v + 72 v 7 ^ (5) 
with integer 71,72- 

Notice that ([5]) is not an identity in the initial vari- 
ables m„ Hi because an integer can have several differ- 
ent presentations as a sum of two squares, for instance, 
4-tuple {{-l,4},{2,-5},{-4,l},{5,-2}} is an example of II- 
type solution, with all weights ji = 1 and g-class indexes 
qi = 17,^2 = 29. 

The I- and H-type of solutions describe substantially dif- 
ferent energy exchanges in the fc-space. The II-type res- 
onances are called angle-resonances [23I ] and consist of 
wavevectors with pairwise equal lengths, i.e. |ki| = | IC3 1 
and |k 2 | = |k 4 | or |ki| - |k 4 | and |k 2 | - |k 3 |. Thus, 
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these resonances do not transfer energy outside of the 
initial range of |k| and, therefore, cannot provide an en- 
ergy cascade mechanism. However, these resonances can 
redistribute energy among the initial wavenumbers, in 
both the direction k/|k| and the scale |k|. Since the ini- 
tial support of energy in |k| cannot change, the II- type 
resonances alone would form a finite dimensional system, 
and it would be reasonable expect a relaxation of such a 
system to a thermodynamic Rayleigh- Jeans distribution 
determined by the initial values of the motion integrals 
(the energy and the waveaction in the this case). Note 
that such a thermalization could happen only among the 
resonant wavenumbers, and many modes which are ini- 
tially excited but not in resonance would not evolve at 
all. (Such an absence of the evolution was called " frozen 
turbulence" in Q where the capillary waves were stud- 
ied for which there are no exact resonances). Whether 
the thermalization does occur in finite clusters and under 
what conditions (e.g. the cluster size etc.) are interesting 
questions that remain to be studied in future. 

On the other hand, the I-type resonances are called scale- 
resonances [23d . and they can generate new wavelengths, 
- this corresponds to 3 or 4 different weights 7 in Q 
[23j ]. Thus, they are the only kind of resonances that can 
transfer energy outside of the range of initial |k|. 

Before studying the structure of resonances let us notice 
the following simple but important fact. Suppose a quar- 
tet 



{ki,k 2 ,k 3 ,k 4 } 



(G) 



is a solution of ([3]), then each permutation of indexes 
I <-> 2, 3 <-> 4 or simultaneous 1 <-> 3 and 2 <-> 3 will gen- 
erate a new solution of ((6]); in general case, all together 
8 different symmetry generated solutions. Of course, in 
some particular cases, when some of the vectors belong- 
ing to a quartet coincide, the overall number of symmetry 
generated solutions can be smaller. For instance, the 
quartet {{{0, -54}, {0, 294}, {90, 120}, {-90, 120}} 
is part of the cluster of 8 symmetry 
generated solutions, while the quartet 
{{17, 31}, {-153, -279}, {-68, -124}, {-68, -124}} 
belongs to the cluster of 4 such solutions. On the Fig. [T] 
the smallest tridents cluster is shown, with and without 
multi-edges. Graphical presentation of a cluster as a 
graph with multi-edges is of course mathematically 
correct but would make the pictures of bigger clusters 
somewhat nebulous (see Fig. [H left panel) . Notice that 
although formally we have solutions with multiplicities 
due to the symmetries, physically all these solutions 
correspond to the same quartet, and therefore we count 
them as one. Therefore, further on we omit multi-edges 
in our graphical presentations as it is shown in Fig. [1] 
right panel. 

In the next sections, on all the Figures, the structure 
of resonance clusters in the spectral space is shown as 





Figure 1: Color on line. The smallest tridents cluster formed 
by 8 symmetry generated solutions, two presentation forms 
are given: with multi-edges (left panel) and without multi- 
edges (right panel) 



follows: each wave vector is presented as a node of integer 
lattice and nodes belonging to one solution are connected 
by lines. 



SCALE- RESONANCES 



We have studied the cluster structure of the scale- 
resonances in the spectral domain |m|, |n| < 1000. In our 
computational domain we have found 230464 such reso- 
nances, among them 213760 collinear (i.e. all 4 wave vec- 
tors are collinear) and only 16704 (7.25%) non-collinear 
resonances. However, the nonlinear interaction coeffi- 
cient in collinear quartets is equal to zero and, there- 
fore, 



they have no dynamical significance [24|. On the 
other hand, for mathematical completeness we will con- 
sider these solutions too, because there may exist 4- wave 
systems with the same dispersion law but with different 
forms of interaction coefficients such that are not neces- 
sarily zero on the collinear quartets. 




Figure 2: Color online. Structure of collinear resonances in 
the spectral domain |k| < 100 (left panel) and of non-collinear 
resonances in the same domain (right panel) 

In the Table I the structure of all clusters is presented 
while in the Table II - the structure of non-collinear is 
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given; cluster length is the number of quartets belonging 
to one cluster. 



Cluster 


Number of 


Cluster 


Number of 


length 


clusters 


length 


clusters 


1 


43136 


56 


3 


2 


1256 


60 


6 


8 


452 


72 


1 


10 


184 


80 


1 


12 


20 


92 


2 


16 


14 


128 


1 


20 


10 


152 


1 


40 


7 


176 


1 


48 


2 


208 


1 



Table I: Clustering in the entire set of quartets in the domain 
1000x1000 (symmetrical solutions are omitted) 



Cluster 


Number of 


Cluster 


Number of 


length 


clusters 


length 


clusters 


1 


1312 


10 


18 


2 


48 


12 


6 


3 


8 


16 


8 


6 


14 


34 


2 


8 


4 


46 


2 



Table II: Clustering in the non-collinear subset of quartets in 
the domain 1000x1000 (symmetrical solutions are omitted) 

In Fig. [^structure of collincar and non-collinear quartets 
is shown ia a smaller domain y mf + n 2 < 100. 



Collinear quartets 

First of all, let us make an important remark. If a 
4-tuple {{toi, ni\, {m 2 , fi2}, {m 3 , n 3 }, {7714, ^4}} consists 
of all collinear wave vectors, then the ratio |mi|/|?ii| is 
the same for all 4 wave vectors. Let us assume that 
rrii, rii ^= 0, then 

^ k|/K| = c, Vz = 1,2,3,4, (7) 

where c is an arbitrary finite rational constant and 

|k l | 1 / 2 = (TO?+ ^ ) l/4 =TO l/2 (1+c2) l/4 ] (8) 

and Sys.Q takes the form 

iW^ + ma 1 / 2 = m 3 1 / 2 + m 4 1 / 2 , 
I mi + m 2 = m 3 + , 

[(m 1 m 2 ) 1/2 = (m 3 TO 4 ) 1 / 2 =*> |mi| = |m 3 m 4 |/|m 2 |. 

Now we can compute mi taking arbitrary integer 
m 2 , m^i provided that |m 3 ?7i4| is divisible on \m2\1 and 
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Figure 3: Color online. First non-axial trident 

keeping in mind that n, = cffij, we can find all collincar 
solutions with c 7^ 0. Obviously, a rational number c de- 
fines a line in the spectral space and not all the lines 
are allowed. Case c = corresponds to the solutions 
lying on the axes X (n^ = crrii) and Y (m, = crii), i.e. 
with all Tfii — or n% — correspondingly, e.g. 4-tuple 
{{4, 0}, {-49, 0}, {-36, 0}, {-9, 0}}. Parametrization of 
the resonances in this case was first given in 24J. As 
we have already mentioned before, these quartets are dy- 
namically irrelevant because there is no nonlinear inter- 
action within these quartets [24| . 



Non-collinear quartets 



Tridents 



Non-collinear scale-resonances have been studied first in 
0, [lU in the spectral domain yj in 2 + nf < 100 and a 
special type of quartets named tridents have been singled 
out. By definition, a wave quartet is called a trident if: 

1) There exist two vectors among four in a quartet, say 
ki and k 2 , such that ki fj k 2 , thus they satisfy 

(ki • k 2 ) = -fcifc 2 ; 

2) Two other vectors in the quartet, k 3 and k4, have the 
same length: A: 3 = ; 

3) k 3 and k4 are equally inclined to ki, thus 
(^ • k 3 ) = (ki ■ k 4 ) . 

The followingpresentation for a trident quartet has been 
suggested in [3, HU , 

ki = (a, 0), k 2 = (-&, 0), k 3 = (c, d), k 4 = (c, -d) , (9) 

and two-parametric series of solutions has been written 
out: 



(s 2 +t 2 + st) 2 , b={s 2 +t 2 -st) 2 
2st{s 2 + t 2 ), d= s 4 -t 4 , 



(10) 



with arbitrary integer s, t. It is easy to check that vectors 
ki, k 2 , k 3 , k4 belong to the same class Cl\, with weights 



7i 



s 2 + t 2 + st, 



",2 



s 2 +t 2 - st, 7 3 



74 



s 2 + t 2 , 
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and obviously 71 + 72 = 73 + 74, Vs, t6Z. 

Parameterization ([9]) corresponds to the tridents ori- 
ented along the X-axis with its vectors ki and kz and, 
therefore, we will call them axial tridents. There ex- 
ist also non-axial tridents, for instance, the quartet 
49, 49}, {-9, -9}, {5,35}, {35,5}}. As mentioned in 
,115(, all the non-axial tridents can be obtained from 
the axial ones, ([9]), via a rotation by angles with rational 
values of cosine combined with respective re-scaling (to 
obtain an integer-valued solution out of rational-valued 
ones). 



Cluster length 


Number of clusters 


1 


1320 


2 


48 


6 


40 


12 


10 


18 


2 


22 


2 



Table III: Clustering in the subset of tridents in the domain 
1000x1000 (symmetrical solutions are omitted) 

In the computational domain |m|,|n| < 1000, we have 
found 13888 non-axial tridents, the first non-axial tri- 
dent is shown on the Fig. O Among 13888, 13504 
tridents have no vectors on any axis and 384 tri- 
dents have just one vector on an axis (for instance, 
{{180, 135}, {0,64}, {120, 119}, {60, 80}}). The total 
amount of all possible tridents is 14848, thus only 960 
are axial (6.5% of the total number). The data on tri- 
dents' clustering are given in the Table III. 
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Figure 4: Color online. First non-trident 
Clusters 

In the previous Section we have shown that there are 
three different types of scale-resonances: collincar quar- 
tets, tridents and non-tridents. There are also clus- 
ters formed by different types of scale-resonances, for 
instance, clusters containing tridents and non-tridents 
(Fig. O left panel) or collinear and non-collinear quar- 
tets (Fig. [5l right panel). Notice that, since collinear 
quartets of gravity water waves have zero interaction co- 
efficients [24| . the cluster shown on the right panel can 
dynamically be regarded as two independent quartets. 
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Non-tridents 



Our study of the resonance solution set in the spec- 
tral domain \rrii\, \ni\ < 1000 shows that not all 
non-collinear cascading quartets are tridents. They 
are called further non-tridents, e.g. a quartet 
{{990, 180}, {128, 256}, {718, 236}, {400, 200}} is a non- 
trident quartet. 

The overall number of tridents is 14848 while the number 
of non-tridents is 1856. Notice that the first non-trident 
quartet {{180, 135}, {0, 64}, {120, 119}, {60, 80}} lies in 
the spectral domain |k| < 225. This means, that if we 
are interested only in the large-scale quartets, say, quar- 
tets with |k| < 100, the complete set of scale-resonances 
consists of 1728 quartets, among them - 1632 collinear 
quartets and 96 tridents, but no non-tridents yet. 



Figure 5: Color on line. Left panel: A shortest cluster 
formed by both tridents and non-tridents; cluster length is 
8 (64 with 8 symmetries), among them 6 (48 with 8 sym- 
metries) tridents and 2 (16 with symmetries) non-tridents. 
Right panel: A shortest cluster formed by both collinear 
and non-collinear quartets; cluster length is 8 (48 with 8 and 
4 symmetries), among them 6 (32 with 8 and 4 symmetries ) 
collinear and 2 (16 with 8 symmetries) non-collinear quartets. 
No multi-edges are shown. 

One of the most important characteristics of the res- 
onance structure is the wavevector multiplicity (intro- 
duced in (20l|). which describes how many times a given 
wavevector is a part of some solution. In the Table IV 
the wavevectors multiplicities are given for non-collinear 
quartets. It turned out that 91% of all these wavevec- 
tors (6720 from overall amount 7384) have multiplicity 
1 (counting the 8 symmetry generated solutions as the 
same quartet.) 
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Multiplicity 


Amount of vectors 


1 


6720 


2 


424 


3 


192 


4 


36 


5 


8 


6 


4 



Table IV: Wavevectors multiplicity computed for the non- 
collinear quartets in the spectral domain |rrij|,|7ii| < 1000 
(symmetrical solutions are omitted) 



In Fig. [TJ we see an example of such a simple cluster con- 
sisting from just one physical quartet. In the Fig. [5] the 
cluster of 4 connected quartets is shown (the symmetry 
generated solutions are omitted), with all together only 
7 different wave-frequencies. 




Figure 6: Color on line. Example of a non-trident cluster of 
length 4. 



Obviously, the cluster structure defines the form of dy- 
namical system corresponding to the cluster. The main 
motivation of our detailed study of clusters is, of course, 
in constructing an isomorphism (i.e. one-to-one corre- 
spondence) between a cluster and a dynamical system. In 
[2l| this construction has been presented for an arbitrary 
3-wave resonance system, with triads as primary elements 
of the planar grap h. Its implementation in Mathematica 
was given in 



221 ] . where interaction coefficients similar 
to Z in (jTHJ) were also computed. To construct this iso- 
morphism two following facts were used: 1) in a 3-wavc 
resonance system only scale-resonances exist, and 2) if we 
add an arbitrary triad to a cluster, in general we always 
add some new wave- freq uencies as well (the only excep- 
tion is identified in [21(). A 4- wave resonance system 
does not possess these nice properties; on the contrary - 
there exist scale- and angle-resonances, and most of the 
quartets are parts of symmetry generated solution sets. 
It would be a challenge to develop a general approach to 
construct dynamical systems for resonance clusters in an 
arbitrary 4-wavc system. 



ANGLE- RESONANCES 



Until now we have studied the structure of the scale- 
resonances which are described by ((4|) and which gen- 
erate new scales (i.e. new values of |fc|, thereby provid- 
ing the energy cascade mechanism. The angle-resonances 
represented by ((5|) do not generate new scales but they 
can redistribute energy among the modes which were ex- 
cited initially (or which are forced externally). Taken 
on their own, these resonances could lead to the ther- 
mal equilibrium distribution of the energy and the wave- 
action among the (initially excited) resonant waves, if 
the number of such waves is large, or they could lead 
to a periodical behavior or a strange attractor, if the 
number of the initially excited modes is small. The im- 
portant fact, however, is that the scale- and the angle- 
resonances are not independent and can form a mixed 
cluster containing both types of this resonances [23] ]. 
The energy cascade mechanism in such a mixed cluster 
is presented schematically in FigJT] where the quadran- 
gles Si and S2 denote scale-resonances so that 4- 
tuples (Vi t i, Vi,4) and (V^i, V2, 4) represent scale- 
resonances and squares A\, ...Ai, .., A n represent angle- 
resonances. 




. n,3 [ 

.0. 



Figure 7: Color on line. Schematic presentation of mixed 
resonance cascade 



If a wave takes part simultaneously in angle- and scale- 
resonances, corresponding nodes of S±,S2 and Aj are 
connected by (red) dashed arrows so that V\ t \ = V n .2 
and V nt i = V<!,3- We have found many examples of such 
mixed clusters in our solution set, for instance: 

Si = {{-64, -16}, {784, 196}, {144, 36}, {576, 144}}, 

A n = {{-64, -16}, {4, 16}, {-64, 16}, {4,-16)}}, 
with V h i = V n ,2 = (-64, -16), 

5*2 = {{-49, -196}, {4, 16}, {-36, -144}, {-9, -36}}, 
which is further on connected with an angle-resonance 
A h = {{-49, -196}, {784, 196}, {-49, 196}, {784, -196}}, 
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(not shown in FigfTJ) via a node-vector (—49, —196). 

The energy flux over scales due to the mixed cluster is 
weak: one wave can participate in a few dozen of angle- 
resonances (see Fig. [5]), which means that only a small 
part of its energy will go to a scale-resonance. More- 
over, even this weak energy cascade may terminate at 
finite wavenumber if it happens to move along a finite 
cluster which terminates before reaching the dissipation 
range (or before reaching the range of high frequencies 
where the nonlinearity gets large enough for the quasi- 
resonanccs to take over the energy flux from the exact 
resonances). 




Figure 8: Color on line. Wave (64,0) takes part in 2 scale- 
resonances, both non-tndents. The upper number in the circle 
is m, the lower is n, and (red) thick lines drawn between 
vectors on the same side of the Eqs.(|3}. Wave (119,120) takes 
part in 1 scale-resonance and in 12 angle-resonances. 

The number of angle-resonances in the spectral domain 
\ni\ < 1000 is of order of 6 • 10 8 while the number of 
scale- resonances in the same domain is of order of 8 • 10 5 , 
among them less than 2 • 10 4 are non-collinear and do play 
a role in the energy exchange among the modes within 
the quartets (see next the Section). It is easy to see 
that an arbitrary wavevector (m, n) takes part in infinite 
number of resonances if spectral domain is unbounded. 
Indeed, let us fix to and n, then a quartet 



(m, n)(t, —n) — ► (to, —n)(t, n) 



(11) 



a scale-resonance with arbitrary t = 0, ±1, ±2, . . .. More 
involved 5-parametric series of angle-resonances found 
from the following considerations. For angle-resonances 
of four wavevectors (a, 6)(c, d) — ► (p, q)(l, m), Eqs.lJHJ can 
be rewritten as 



f b 2 = p 2 + q 2 7 
c = p + l, b 



c 2 + d 2 = I 2 
d = q + to 



(12) 



Simple algebraic transformations and known parameteri- 
zations of the sum of two integer squares (e.g. for a circle 
or for the Pythagorean triples) yield 



t 2 )/(s 2 +t 2 ), 
2^ K tl i „2\ 



= (s 2 

P=(/ 2 -<? 2 )/(/ 2 + <r 

d = (a 2 +b 2 + ac — ap 



b = 2st/(s 2 + t 2 ), 

q = 2fg/(f 2 + g 2 ), (13) 
cp-bq)/(q - b) 



This is an easy task to check then that the solutions of 
(|12p can be written out (perhaps with repetitions) via five 
integer parameters s, t, /, g, c (rational solutions should 
be renormalized to integer). Notice that (fTTj) degenerates 
to trivial resonances 



(m,0)(t,0)- (m,0)(t,0), 



(14) 



if n = 0, i.e. it does not include any resonances of 
wavevectors of the form (to,0), for instance (1,0). In 
this case the choice / = 1, g = 1 in (Ti"3|) gives 



(l,0)(c,l + C )^(0,l)(l + c,c). 



(15) 



Analytical series are very helpful not only for computing 
resonance quartets and clusters structure but also while 
investigating the asymptotic behavior of interaction co- 
efficients. 




10 20 



40 50 60 70 



Figure 9: The multiplicities histogram for angle-resonances. 

The multiplicity histogram for the angle-resonances is 
shown in Fig[9j On the axis X the multiplicity of a vec- 
tor is shown and on the axis Y the number of vectors 
with a given multiplicity. This graph has been cut off - 
multiplicities go very high, indeed the vector (1000,1000) 
takes part in 11075 solutions. All this indicates that the 
angle-resonances play an important role in the overall 
dynamics of the wave field. 



DYNAMICS 



The wave field 



Once the clusters are found, one can consider an evolu- 
tion of amplitudes of waves that belong to each individual 
cluster by considering a respective reduction of the dy- 
namical equation. For the gravity wave case, the appro- 
priate dynamical equation is Zakharov equation for the 
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complex amplitude a^it) corresponding to the k-modc, 
.dak 



dt 



E 

k1.k2.k3 



nk,ki 
k 2 ,k 3 l 



Tu z C- a kl ak 2 ak 3 e 



(16) 

where T ka k = T(k, ki, It2, k 3 ) is an interaction coeffi- 
cient for gravity water waves which can be found in [25| . 

For very weak waves, amplitudes Ok(t) vary in time 
much slower than the linear oscillations. The factor 
e *( Wk + Wkl - Wk2 - Wk3 )t Qn RHg wiU rap idly oscillate for 

most waves except for those in an exact resonance for 
which this factor is 1. Thus, only the resonant modes 
will give a contribution to the dynamics in this case and 
the oscillating contributions of the non-resonant terms 
will average out in time to zero and will not give any 
contribution to the cumulative change of cik(t). Leaving 
only the resonant terms, we have 



l ~dT 



R 

^.ka^ak^ks, 

k1.k2.k3 



(17) 



where ^ means a summation only over kx , k 2 and k 3 
which are in resonance with k. Obviously, we should 
consider k's from the same cluster only (i.e. solve the 
problem for one cluster at a time). Note that the fast 
timescale of the linear dynamics completely disappeared 
from this equation. Thus, paradoxically, the dynamics of 
very weak waves in finite boxes is strongly nonlinear: it 
is more nonlinear than in WTT which works for larger 
amplitudes and where quasi-resonances ensure phase ran- 
domness. This explains the fact found in the three-wave 
example that even relatively large clusters often exhibit 
a periodic or quasi-periodic behavior (26[. 

A study of the system (|17p is possible analytically for 
small clusters (perhaps even integrating the system in 
some lucky cases) and numerically for large clusters, 
which is an interesting subject for future research. Some 
properties, however, can already be seen in the example 
of a single quartet which has been studied before 27 1. 
Let us briefly discuss these properties. 



A quartet 

The dynamical system describin g sl owly changing ampli- 
tudes of a quartet has the form [271 ] : 

idi = iZa^a^a^ 
ici2 = 2Za\a^a^ 
iaj, = 2Z* a\aia,2 
104 = 2Z* 0,30,10,2 



(18) 



nkl.k2 

k 3 ,k 4 ' 



The mathematical analysis of the sys- 
tem ()18[) can be performed similar to what has been done 



in [28[ for an intcgrablc 3-wave system of resonantly in- 
teracting planetary waves, though computations of the 
modulus of elliptic integral are more involved and variety 
of different dynamical scenarios is substantially readier 
(see 2JJ for details). The general answer can be given in 
terms of Jacobean elliptic functions. 

Obviously, the quartets with interaction coefficient Z = 
do not influence the general dynamics of the wave field 
at the corresponding time scale. As it was shown in [2~1 |. 
for all collinear quartets Z = 0, and therefore they can 
be excluded form consideration. 

On the other hand, in [29[ so-called "degenerate quar- 
tets" (tridents in our terminology) have been studied and 
it was established numerically that they have strictly pe- 
riodic behavior: after appropriate translation in the hor- 
izontal plane, two snapshots of the free-surface taken at 
t = and t = T are identical. 

For a general quartet and for arbitrary and arbitrary ini- 
tial conditions, the system (fl8|) "does not exhibit strict 
periodicity" in numerical simulations (27| . Thus, the gen- 
eral formulae for the solutions of (|18|1 have to be stud- 
ied in more details in order to distinguish periodic and 
non-periodic dynamics of an arbitrary quartet. On the 
other hand, since our main interest in this paper is the 
large-scale dynamics, all scale-resonances with a non-zero 
interaction coefficient are tridents and therefore demon- 
strate a periodic time behavior. Of course, if a trident is 
involved into a cluster with some other quartets then one 
should be cautious about the predictions obtained for an 
isolated quartet. 



SUMMARY AND DISCUSSION OF RESULTS 

• In this paper, we studied properties of deep water grav- 
ity waves bounded by a square periodic box. At very 
small wave amplitudes, when the nonlinear resonance 
broadening is less than the fc-space spacing, WTT fails 
and only the waves which are in exact four-wave reso- 
nance can interact. This situations appears to be typical 
for all existing numerical simulations [9( and laboratory 
experiments [lOj. Thus, to understand the wave behavior 
in laboratory experiments and in numerical simulations 
it is crucial to study exact resonances among discrete 
wave modes, which was the focus of the present paper. 
Of course, the notation of "very small amplitudes" has 
to be worked out explicitly for interpreting the results of 
laboratory experiments. The smallness of amplitudes is 
defined by the choice of a small parameter < e <C 1 
and depends of the intrinsic characteristics of the wave 
system, for instance, for atmospheric planetary waves it 
is usually taken as the ratio of the particle velocity to 
the phase velocity which allows to obtain explicit esti- 
mation [30I ] for a wave amplitude a(m,n) (correspond- 
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ing to the weakly nonlinear regime) as a function of to 
and n. For the water surface waves, the wave steepness, 
\a(m, n)\/ (to 2 + n 2 ) 1 / 2 ^, is usually taken as a small pa- 
rameter e and e ~ 0.1 corresponds then to the weakly 
nonlinear regime. For such weakly nonlinear waves to 
feel discreteness of the fc-space, their nonlinear frequency 
broadening has to be less than the distance between ad- 
jacent fc-modes. In terms of the wave s teep ness this con- 
dition reads e < (m 2 + n 2 ) -1 / 8 , see 



e stecpr 



• Wc found numerically all resonant quartets on the set of 
1000 x 1000 modes making use of the g-class method orig- 
inally developed in [l8[ . We found that all resonant quar- 
tets separate into noninteracting with each other clus- 
ters. Each cluster may consist of two types of quartets: 
scale- and angle-resonances. The angle-resonances can- 
not transfer energy to any /c-modes which are not already 
present in the system. They cannot carry an energy flux 
through scales and their main role is to thermalize the 
initially excited modes. The scale-resonances are much 
more rare than the angle ones and yet their role is im- 
portant because they are the only resonances that can 
transfer energy between different scales. Most of the 
scale-resonances, but not all, are of the trident type, for 
which a partial parametrization can be written out ex- 
plicitly. If one is interested in large-scale modes only, say 
100 x 100 domain, then the tridents are the only scale- 
resonances, which is very fortunate due to the available 
parametrization. 

• Even though the angle-resonances cannot cascade en- 
ergy, they are important for the overall cascade process 
because they are involved in the same wave clusters with 
the scale-resonances. One wave mode may typically par- 
ticipate in many angle resonances and only one scale res- 
onance. Thus, one can split large clusters into "reser- 
voirs" , each formed by a large number of angle quartets 
in quasi-thermal equilibrium, and which are connected 
with each other by sparse links formed by scale quar- 
tets. This structure suggests a significant energy cascade 
slowdown and anisotropy with respect to the infinite-box 
limit. Further study is needed to examine the structure 
of such large clusters and possible energy cascade routes 
from the region of excitation at low wavenumbers to the 
dissipative large k range. In the essentially finite domain 
the situation is opposite. Quite recently results of the 
laboratory experiments with surface waves on deep water 
were reported [31 . in which regular, nearly perma- 
nent patterns of the water surface have been observed. A 
feasible way to interpret these results would be 1) to es- 
tablish that the conditions of the experiments correspond 
to the weakly nonlinear regime; if yes - to proceed as 
follows: 2) to compute all exact resonances in the wave- 
lengths range corresponding to those in the experiments; 
3) to demonstrate that for chosen wave-lengths and the 
size of laboratory tank no scale-resonances appear, 4) 
to attribute the regular patterns to the corresponding 



angle-resonances. Obviously, the scale-resonances would 
produce the spectrum anisotropy and disturb the regular 
patterns. 



• We also discussed consequences of the cluster structures 
for the dynamics, and argued that one should expect a 
less random and more regular behavior in the case of 
very low amplitude waves with respect to larger (but still 
weak) waves described by WTT. More study is needed 
in future both analytically, for small clusters, and nu- 
merically, for large clusters. A particularly interesting 
question to answer in this case is about any possible uni- 
versal mechanisms of transition between the regular dy- 
namics to chaos and possible coexistence of the regular 
and chaotic motions. 



• Last not least. The knowledge of the 2D-resonancc 
structure might yield new insights into the origin of some 
well-known physical phenomena, for instance, Bcnjamin- 
Fcir (B-F)instability [H| or McLean instability. This is 
"a modulational instability in which a uniform train of 
oscillatory waves of finite amplitude losses energy to a 
small perturbation of waves with nearly the same fre- 
quency and direction" [35j . As it was shown recently 
H, 37 1, the modulational instability, though well 



established not only water waves theory but also in plas- 
mas and optics, has to be seriously reconsidered. It 
turned out that 1) it can be shown analytically that arbi- 
trary small dissipation stabilizes the B-F instability, and 
2) results of laboratory experiments show that B-F the- 
ory generally over-predicts the growth rate. Moreover, 
the growth rate changes with the time j3{|. Some re- 
searchers state even that "... this effect is far less signif- 
icant than was believed and should be disregarded" [38( . 
The other way to treat the problem would be to try and 
and explain the modulational instability through non- 
collinear (that is, essentially two-dimensional) exact res- 
onances [39j. Similar questions arise in the study of 
McLean instabilities defined by the magnitudes of the 
water depth on which surface waves are studied. For in- 
stance, as it was demonstrated in [4pj ] . in some regimes 
of shallow water the instabilities are due to higher order 
resonances among 5 to 8 waves. It would, therefore, be 
interesting to see how the B-F and McLean instabilities 
are modified by the finite flume effects and the corre- 
sponding discreteness of the wave resonances, and to see 
what role this could have played in the past laboratory 
and numerical experiments. Results presented in our pa- 
per can be regarded as a necessary first step for such 
an investigation, and the future work would involve ap- 
plication of the g-class method to computing the higher 
order resonances which may be involved in modulational 
instabilities. 
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